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The flow properties of confined vortex matter driven through disordered mesoscopic channels 
are investigated by mode locking (ML) experiments. The observed ML efi'ects allow to trace the 
evolution of both the structure and the number of confined rows and their match to the channel 
width as function of magnetic field. From a detailed analysis of the ML behavior for the case 
of 3-rows we obtain (?) the pinning frequency fp, (ii) the onset frequency fc for ML (oc ordering 
velocity) and (Hi) the fraction Lml/L of coherently moving 3-row regions in the channel. The 
field dependence of these quantities shows that, at matching, where Lml is maximum, the pinning 
strength is small and the ordering velocity is low, while at mismatch, where Lml is small, both the 
pinning force and the ordering velocity are enhanced. Further, we find that fc cc fp, consistent with 
the dynamic ordering theory of Koshelev and Vinokur. The microscopic nature of the fiow and the 
ordering phenomena will also be discussed. 



I. INTRODUCTION 

Vortex arrays (VA's) in type II superconductors are 
exemplary systems to study non-equilibrium states of 
driven periodic media in various pinning environments. 
A particularly interesting phenomenon in this context is 
that of a dynamic transition from an elastic, coherent 
flow state at large velocities to a plastic, incoherent flow 
state at small velocities. The first theoretical descrip- 
tion of this issue was provided by Koshelev and Vinokur 
(KV) [1] and refined in subsequent studies [2-5] which 
predicted various novel flow states, including a moving 
glass characterized by elastically coupled chains oriented 
along the flow direction and a moving transverse smec- 
tic with decoupled flow chains. Such structures and the 
dynamic transitions between them have been extensively 
studied in a number of numerical simulations [5-7]. 

Experimentally, a diversity of flow states has been re- 
ported in direct imaging experiments on NbSe2 crystals 
[8-10]. However, quantitatively the effect of pinning 
strength and/or temperature on the ordering velocity 
has been studied only through dc transport experiments 
[11-13], based on the assumption that an inflection point 
in dc current-voltage {IV) curve (i.e. a peak in differen- 
tial resistance) marks the dynamic transition. Recently, 
different explanations have been given as the reason for 
such inffection point, like macroscopic coexistence of two 
phases [14] and a change in the self organized, large scale 
morphology of vortex rivers [15]. Thus, a more direct, 
microscopic probe is required to study systematically the 
velocity, magnetic field and temperature dependence of 
dynamic ordering. 

Recently, we reported on the use of mode-locking (ML) 
experiments as a direct probe of ordering [16]. The ML 
phenomenon occurs due to coupling between, on the one 
hand, collective lattice modes of frequency fmt — qvdc/a, 



with q an integer and a the lattice periodicity, which 
occur when a VA moves coherently with velocity Vdc 
through a pinning potential [10,17], and, on the other 
hand, a superimposed rf-drive of frequency / at an inte- 
ger fraction 1/p of fint- This coupling produces steps 
in the dc-transport {IV) curves when Vdc = {p/q)!^, 
[7,18-21], similar to ML steps in sHding charge density 
waves (CDW's) [22,23] and giant Shapiro steps in Joseph- 
son junction arrays [24]. However, on decreasing the ve- 
locity Vdc or increasing the temperature, incoherent ve- 
locity fluctuations due to quenched and thermal disorder 
suppress the collective lattice mode and reduce the width 
of the ML steps. Approaching the regime of fully plas- 
tic or liquid flow, the ML amplitude eventually vanishes 
[7,16] and the ML frequency fc at which this occurs pro- 
vides a direct measure of the ordering velocity: Vc = fcO-- 

The particular system of our studies consists of meso- 
scopic flow channels in a disordered, strong pinning envi- 
ronment [21,25]. The geometry of the samples is sketched 
in Figs. l(a),(b). Vortices inside the channels are con- 
fined by strongly pinned vortices in the channel edges 
(CE's). When a force is applied along the channel, 
the shear interaction with these CE vortices provides 
the dominant pinning mechanism impeding the channel 
fiow. Further, since the natural lattice (row) spacing is 
Go — 1.075-\/<i>o/i? {bo — V^ao/2), on varying the mag- 
netic field B one can go through a series of structural 
transitions from n to nil vortex rows in the chan- 
nel (typically n < 10). Besides its relevance for the 
study of (dynamic) structural transitions of vortex mat- 
ter in quenched disorder, the physics of this system is also 
closely related to layering transitions in confined fluids, 
fiow of colloids in mesopores and mesoscopic friction. 

In Ref. [21] we have given a short account of the in- 
teresting phenomena which this system displays. Firstly, 
due to the structural transitions, the dc-critical current 
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for channel flow (yield strength) oscillates with field as 
shown in Fig. 1(e). At a given field, the dc IV (force- 
velocity) curve with superimposed rf-current exhibits the 
ML eflfect, as shown in Fig. 1(c). The ML condition 
in this case attains a form which is particularly useful 
to study the structural transitions. The voltage at 
which the fundamental(p = q = 1) ML step occurs is 
given by [21]: 



(1.1) 



with $0 the flux quantum, Nch the number of channels 
measured simultaneously and n the number of coherently 
moving chains in each channel. 
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FIG. 1. (a) The channel device (side view) consisting of 
strong pinning NbN (dark) and weak pinning amorphous 
NbGc. The current and ficki directions arc indicated, (b) 
Sketch of the vortex structure around a channel (top view). 
The effective channel width w is also indicated, (c) Typi- 
cal ML step in a dc-IV curve at 70mT and a superimposed 
rf-current of 3MHz. (d) Normalized ML voltage Vi,i/(/$o) 
versus magnetic field, (e) Critical current Ic, determined from 
a 1/xV criterion, versus field. 



plateau, where an n row configuration matches the chan- 
nel width, it is minimum. As shown in Rcf. [21], this phe- 
nomenon is caused by positional disorder (roughness) of 
the vortex configuration in the CE's: first of all, on mov- 
ing away from the matching field this disorder enhances 
transverse fluctuations of vortex chains in the array, im- 
peding the flow. Secondly, close to mismatch, part of the 
(moving) n-chain regions within each channel may switch 
to n ± 1 . In between the n and nil regions quasi-static 
fault zones with misaligned dislocations develop where 
the vortex trajectories are jammed. We note that the 
presence of degrees of freedom transverse to the aver- 
age flow velocity, both in our system and vortex lattices 
in general, forms an important difference with CDW's. 
Particularly, for CDW's the 'displacement' (phase) field 
is a scalar (j){x) and the velocity oc dt(t){x) represents lon- 
gitudinal motion only. Our system, when compared to 
regular vortex lattices, has the unique property that the 
transverse response can be tuned. 

In this paper we study in detail the ML step width as 
function of the rf-amplitude and frequency. The exper- 
iments provide important information on the dynamic 
coherence of the arrays and how this coherency varies 
with mismatch and fiow velocity. We focus on ML in the 
regime where n = 2 — > 3 — > 4, but our findings are repre- 
sentative for other transitions with limited n < 10. The 
paper is organized as follows. In Sec. II we first review 
previous theoretical results on ML at high frequency and 
then present simidation results of a one dimensional (ID) 
vortex chain to show the full frequency dependence of ML 
in a system with only elastic deformations. In Sec. Ill 
we describe the details of our sample and the experimen- 
tal procedure. The experimental results are presented in 
Sec. IV. We find clear evidence for an ordering frequency 
fc below which no ML occurs. Furthermore, from the ML 
data we extract the pinning frequency fp and the total 
length Lml of coherently moving 3-row regions. These 
quantities systematically change with magnetic field: at 
matching, where Lml is large, fp and fc are both small, 
while at mismatch where Lml is small, both fp and fc 
are enhanced. We find that fc fp, independent of mag- 
netic field. In Sec. V we compare these results with the 
ordering theory of Koshelev and Vinokur, and discuss the 
implications for the microscopic nature of the flow and 
the ordering phenomena. 



II. THEORETICAL CONSIDERATIONS 



As observed in Fig. 1(d), on changing field, Vi^i in- 
creases as a staircase, directly reflecting the evolution of 
n with field . A comparison with Fig. 1(e) shows that at 
mismatch fields, i.e. where a transition n ^ n± 1 occurs 
and n and nil steps may coexist, the yield strength 
a Ic is maximum, whereas around the center of the Vi^\ 



The velocity (or frequency) dependence of the ML step 
width is particularly useful as a direct probe of dynamic 
ordering. At present, the ML step width has been stud- 
ied theoretically only in the high frequency limit, where 
perturbation theory allows to obtain an analytical de- 
scription [26-28]. 
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A. Equation of motion 

The 2D displacement field u{r, t) of an elastic vortex 
lattice driven through a pinning environment by com- 
bined rf and dc external forces is (at T = 0): 

7^=Fd+Fr + Fp (2.1) 

with 7 = $o-B//0/ the friction coefficient with pf the flux 
flow resistivity. Fd is a driving force per unit length con- 
sisting of dc and rf terms; |Fd|= jdc^o + jrf^ocos{27rft) 
with jdc and jrf the dc and rf current densities, re- 
spectively. Fr is the elastic restoring force given by 
($o/-B)[(cii - C66)V(V • u) + ceeV^u] with cn and cee 
the compression and shear moduli, respectively [1,5]. In 
absence of the pinning force Fp , the lattice is undistorted 
and flows uniformly: du/dt = Fn/j = V4c+VacCOs{2'Kft) 
with an ac velocity Vac = jV/^o/T) i-e- proportional to 
the rf-drive. 



B. Amplitude of the ML-interference step 

At high velocities where the friction term dominates 

the pinning term in Eq.(2.1), one can treat the pinning 
as a perturbation with respect to the undisturbed rf-dc 
velocity |Fd|/7. We distinguish two cases, namely a pe- 
riodic pinning potential and a random pinning potential. 
In case of a periodic potential with periodicities equal 
to that of the lattice, clastic deformations are absent 
(Fr,=0) and the whole lattice behaves as a single particle 
with overdamped dynamics in a sinusoidal potential. At 
large drive this case is described analogous to a voltage 
biased, resistively shunted Josephson junction [29]: sub- 
stituting u = vtin Eq.(2.1) and assuming Fp = sm{ku) 
with /i the maximum slope of the potential and k = 27r/a, 
one can show that as first order correction a step anomaly 
appears in the dc velocity-force characteristics at the ML 
condition Vdc = paf- The current density width of the 
pth step oscillates with the rf-drive according to 

Ajp.i = 2jc\Jp{vac/fa)\ (2.2) 

with jc = /i/*&o the critical current density and Jp the 
Bessel function of the first kind of order p. Note that no 
subharmonic ML (q > 2) occurs in this model. 

Turning to a VL in a purely random potential, the first 
order perturbation correction has zero mean. The second 
order correction is the lowest order of perturbation that 
provides the ML step. Taking into account the lattice 
distortions due to the random pinning within the elastic 
limit, Schmidt and Hauger [26] showed that the ML step 
can appear at both harmonic and subharmonic ML con- 
ditions Vdc = {p/Q)cif and that the associated width of 
the current density step is: 

Aip.g = 2jc{qk)j^{qvac/fa) (2.3) 



jc = Y.jc{qk) (2.4) 

Q 

where jdqk) is the component of the critical current den- 
sity related to the Fourier transform of the random po- 
tential correlator at wave vector qk. Thus, for random 
pinning the ML step width exhibits a squared Bessel- 
function oscillation with the rf drive. In the following we 
omit the subscripts p and q in the ML step width since 
we will discuss only the fundamental ML phenomenon for 
p = q = 1 (where q = 1 will be justified in Sec. IV). We 
note here that in our case a can be a frustrated lattice 
spacing opposed to the natural one in [26] . 

C. Frequency dependence 

It is clear in Eqs.(2.2) and (2.3) that the dependence 
on frequency and rf drive only appear in the argument 
^ = ''^ac/fO' of the Bessel functions, irrespective of the 
type of pinning. Choosing z by varying the rf-amplitudc 
such that Ji(z) is maximum, the characteristic maximum 
value of the fundamental ML width at high frequency is 

^jmax,P = l-16ic (2.5) 

Aw,fl = 0.67ie(fc) (2.6) 

for periodic and random pinning, respectively. 

These pcrturbational results are only applicable for fre- 
quencies (much) above the so-called pinning frequency 
fp. For a sinusoidal pinning potential, fp is analogous to 
the characteristic frequency of an overdamped Josephson 
junction and it is given by: 

fp = jc^o/ja. (2.7) 

For random pinning, the fundamental ML step involves 
only the k = 27r/a Fourier component of the pinning 
force, Fp{k) = jc{k)^o since it is responsible for the dy- 
namic lattice mode excited at the washboard frequency 
{fint = Vdc/ a). In this case one can define the pinning 
frequency fp [30] such that the strength of the friction 
term 7a/ in Eq.(2.1) equals Fp{k): 

fp = ~jc{k)^oha. (2.8) 

Below fp no analytical result for Ajmax is available, 
not even for sinusoidal pinning [31]. In this regime nu- 
merical simulations are a useful tool to obtain the the- 
oretical value of the fundamental ML width [32,33]. To 
obtain Ajmax vs. / for the case of completely clastic mo- 
tion in our channel system, we have performed molecular 
dynamics simulations of an rf-dc driven ID elastic vor- 
tex chain both in a channel with periodically configured 
static vortices in the CE's and in channels with strongly 
disordered CE vortex arrangements (see the inset of Fig. 
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2(c) and [34.35] for more details). The channel width 
w ^ bo, i.e. the (average) spacing between the first 
pinned rows is 2bo. Vortex interactions were modelled 
by the London potential with X/uq = 1 (A is the pene- 
tration depth) and the average vortex spacing a in the 
channel was chosen equal to that in the CE's, a = uq. 

For the periodic case, the CE potential is sinusoidal 
and the critical current density jc is given by its max- 
imum slope jc = /i$o- Simulating a chain of limited 
length was sufficient since all vortices behave as a single 
particle. Figure 2 shows an example of the ML step for 
superimposed rf-drive of amplitude Vac/i fa) = 2 and fre- 
quency / ~ 3/p. In Fig. 2(c), we summarize the numer- 
ical results of Ajmax,p versus frequency, represented by 
the solid squares. Here, Ajmax,p is normalized by jc and 
the frequency is normalized by fp defined by Eq.(2.7). 
At high frequency f > fp, Ajmax.p saturates at a fre- 
quency independent value Ajg^p/jc — 1.13 very close to 
the theoretical prediction Eq.(2.5) for periodic pinning. 
For smaller /, Ajmax.p starts to decrease around fp and 
then vanishes linearly with /. The whole frequency de- 
pendence of Ajmax,p is well approximated by an empir- 
ical function 



Aj max — 

Aj,tanh(//0.7/p 



(2.9) 



in which we have omitted the subscript referring to the 
periodic pinning potential. Equation(2.9) is represented 
by the solid line in Fig. 2(c). 

For the disordered ID channel, the CE vortices are 
assigned re latively la rge random shifts over distances 
|d| with -\/((V~dp) = 0.12 with respect to the reg- 
ular lattice configuration [35,36]. We used a chain of 
2000 vortices for which the results have become length- 
independent. Due to the disorder, both the numeri- 
cally obtained threshold force cx j^ and the step width 
Ajmax,R are strongly reduced (by a factor 5, [35]) with 
respect to the ordered case. An example of the simulated 
ML step is shown in Fig. 2(b). There are two distinct 
differences with the periodic case, displayed in Fig. 2(a): 
(<) the sharp corners disappear and (ii) both below and 
above the ML condition the curve in (b) is essentially lin- 
ear with the same slope, with a shift at the ML condition, 
very similar to the experimental result in Fig. 1(c). After 
normalization Aj^axM/j^ and f / fp = laf /{jf^o), we 
plot the simulation results collected for various frequen- 
cies in Fig. 2(c) as the open squares. The saturation 
value Ajmax,R{f > fp) = Ajs^R ~ 0-7jc^ very close to 
the result in the random pinning limit in Eq.(2.6). The 
whole frequency dependence is then again well approxi- 
mated by Eq.(2.9) (dashed line in Fig. 2(c)) in which we 
have now implicitly assumed the subscript 'R' referring 
to the quantities in the random case. 

It is worth mentioning that the results for the disor- 
dered channel were insensitive to small changes in the ra- 
tio a/flQ. Further, we note that both in the periodic and 
the disordered channel simulations ML can be observed 
down to / = 0. We believe this is a direct consequence 



of the fact that vortices in the chain remain elastically 
connected, because other simulations in which also trans- 
verse degrees of freedom and plasticity are allowed, do 
not show this feature [37]. 
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FIG. 2. (a) Normalized force-velocity curves for simula- 
tions of a vortex chain with an rf-drive of frequency / — 3/p 
and amplitude Vac/ {fa) = 2 in a periodic ID-channcl; (b) 
same as (a) in a disordered ID channel, (c) Maximum ML 
current density step width from simulations with rf-drive of 
various frequencies and amplitudes: (■) results for a chan- 
nel with periodic CE vortices; (□) results for disordered ID 
channels. The curves show the empirical relation Eq.(2.9) 
with C = 1.16 (drawn line) and C = 0.69 (dashed Une). The 
inset shows the simulation geometry for a periodic CE con- 
figuration (o) and for a disordered CE configuration (•). 



III. EXPERIMENTAL 

The device consists of a strong pinning layer of poly- 
crystalline NbN film on top of a weak pinning amorphous 
(a-) Nbi_:rGe:r film (x w 0.3). The thickness of the NbN 
and a-NbGe films are 50 and 550nm, respectively. Using 
reactive ion etching with proper masking [38], narrow 
straight channels were etched from the top (NbN) layer 
leaving a 300nm{=dch) thick (a-NbGe) bottom layer, see 
Fig. 1(a). The width and length of each channel are 
230nm and 300/zm(= L), respectively. The spacing be- 
tween adjacent channels is 10 /im. Magnetic field was 
applied perpendicular to the films, inducing a vortex ar- 
ray as schematically shown in Fig. 1(b). The transport 
current was applied perpendicularly to the channel, pro- 
viding a driving force parallel to the channel. 

For the ML measurement, we recorded dc voltage by 
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sweeping the dc current with superimposed rf-current. 
The transmission lines for the rf-current were terminated 
by matching circuits very close to the sample. To avoid 
heating, both the sample and the circuits were immersed 
in superfluid ^He. For consistency, all the data presented 
in this paper were taken after field cooling in which the 
magnetic field was applied above Tc's of a-NbGe (2.68K) 
and NbN(llK) and the sample subsequently cooled to 
T =:1.9K, which is much lower than the vortex lattice 
melting temperature « 2.bK for the fields we studied. 
The ML steps in IV curves are always rounded like in 
Fig. 1(c). For definition of the current step width A/, 
we took the derivative of the IV curve and integrated 
over the ML peak in the differential conductance- voltage 
curve with respect to the flux- flow base line [21]. 



IV. RESULTS 

Our measurements were carried out in the magnetic 
field regime where 3-chain structures exist in the chan- 
nels, ranging from 45 to llOmT, see Fig. 1(d). We 
thus focus on the fundamental ML step characterized by 
Vi,il {^of Nch) = 3 originating from coherently moving 
n = 3 regions. At the borders of our field range coexis- 
tence with n = 2 or n = 4 ML steps may occur. 
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FIG. 3. Current width A7 of the fundamental ML step 
vs. rf current Irf taken at 12MHz and 80mT. The maximum 

value A/maa; IS indicated. Solid and dotted curves display 
AJ oc J'\{z) and A7 oc lJi(2;)|, with z oc Irf, respectively. 

In Fig. 3 we show an example of how the funda- 
mental ML width A/ depends on the rf drive Irf for 
/io-ff = 80mT and / = 12MHz. M{Irf) shows oscilla- 
tory behavior with a maximum value i^Imax in the first 
lobe, as indicated in the figure. A qualitative comparison 
of the data with the theoretical predictions Eq.(2.3) (solid 
line) and Eq.(2.2) (dotted line) shows that the data fol- 
lows more closely a than a | Ji| dependence. At ^Imax 
it is found that Irf = (1.9±0.3)/(ic over a broad frequency 
range between fp and 40 MHz (above this frequency the 
experimental error in 7^/ becomes larger). This value is 



in good agreement with Eq.(2.3) which has a maximum 
at 2 = Vac/ fa = Vachdc ~ Irf/hc = 1-8. It is impor- 
tant to note that while the values of the rf-current might 
appear rather large, the actual vortex displacements due 
to the rf-drive at or below the first maximum in Ji{z), 
are less than 1.8/(27r) ~ 0.3 of the lattice spacing. 

The Jl{z) behavior shows that the pinning potential 
due to the vortices in the CE acts as a random potential 
(RF) [39]. The origin of the RP is the strong positional 
disorder of the vortex lattice in the NbN edge material. 
This disorder has recently been observed in scanning tun- 
nelling microscopy experiments on NbN films [40]. We 
further note that we did not observe subharmonic ML 
steps, i.e. there were no ML steps at Vml = ^f^oNch/Q 
with q > 2. In the context of our RP this seems in 
contradiction to the results of others [18]. But those ex- 
periments have been carried out at relatively low fields 
where the RP has short range correlations on a scale 
^ <C a. Fourier components qk with q>2 are needed to 
describe such short range fluctuations and therefore sub- 
harmonic ML steps are seen. In our case, the RP is due 
to the CE vortices which have average spacing ao — a. 
Consequently the most important Fourier component de- 
scribing this RP is the q = I mode, which explains why 
we do not see subharmonic ML steps. 
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FIG. 4. The maximum current width Almax as a function 
of frequency / at 50mT. The solid curve shows the fit ac- 
cording to the empirical function Eq.(2.9). The dotted line 

shows the linear extrapolation to Almax = for definition 
of the dynamic ordering frequency fc. The inset shows the 
onset behavior of Almax and ordering frequency fc for a field 
of llOmT. The solid curve shows a fit of the high frequency 
data to Eq.(2.9). 

Next we turn to the frequency dependence of Almax- 
Figure 4 shows Almaxif) at 50 mT obtained from mea- 
surements similar to those in Fig. 3 at various frequen- 
cies. As in the numerical results in Fig. 2(c), Al^ax 
saturates at a value A/^ ~ 78 /zA at high frequencies, 
while at low frequencies it decreases monotonically with 
/. A large part of the data is well approximated by 
the empirical function discussed in the previous section. 
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^Imax = ^-fstanh(//0.7/p), and wc can extract the pin- 
ning frequency /p=7.8MHz as the remaining fit param- 
eter. However, at low frequency (i.e. small dc-velocity) 
the data lie significantly below the empirical curve. This 
implies that the vortex motion becomes less coherent due 
to the disordered CE's. On reducing /, Almax vanishes 
almost linearly at a finite frequency fc determined by the 
intersection between the dotted line and the Almax = 
axis. In this regime, the rf-current for which AI exhibits 
its maximum value, starts to saturate at a value ~ Ic 
[41]. The collapse of Almax at finite frequency fc is even 
more clearly visible in the data taken at llOmT shown 
in the inset to Fig. 4. Below fc no ML step appears at 
any rf drive, indicating the complete absence of coherent 
3-row motion. 
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FIG. 5. (a) Vi,i//$o versus field, (b) The pinning fre- 
quency fp (■), the dynamic ordering frequency fc (•), and the 
saturation value of the maximum current width AI, (□) for 
n = 3 as a function of field. The data are obtained from fits to 

the measured AI,nax{f)- (c) The coherently moving fraction 
Lml/L of n = 3 regions vs. field, where Lml is determined 
from Eq.(4.1), with a = ao{B) (□) and a = umBm/B (o), 
where the subscript 'M' refers to the matching field, see text. 
All lines are guides to the eye. 

The above fitting analysis was performed on Al^ax (./) 
data taken at various fields in the n = 3 field regime and 
we extracted fp, A/g and the dynamic ordering frequency 
fc- Wc first discuss the results of fp and Alg as a func- 
tion of magnetic field, shown in Fig. 5(b). As observed, 
fp has a minimum at 70 mT, somewhat below the mid- 



dle of the plateau in Vi^i in Fig. 5(a), and it increases 
on approaching either end of the plateau. Hence, the as- 
sociated pinning current density in the coherent 3-chain 
regions jc{k) is small at 70 mT and increases away from 
70 mT. In fact, the value of jc{k) as determined from 
fp using Eq.(2.8) agrees within 30% with jc as deter- 
mined directly from dc-/F curves. Fig. 1(e). Meanwhile, 
A/g exhibits a field dependence which differs consider- 
ably from that of fp: it has a broad peak around _B = 50 
mT and then decreases with increasing field. Clearly, 
this behavior can not be explained by simply assuming 
Ah{B)(xjc{B)^fp{B). 

At this point we note that theory assumes all vortices 
in the channel are moving coherently. However, in our 
experiment only a fraction of the vortices move coher- 
ently. Specifically, an n-row region may locally break up 
due to the strong edge disorder or it may coexists with 
n ± 1-row regions (with different ML voltages) due to 
mismatch [21]. We define the total length of mode- locked 
regions with 3 coherently moving rows as Lml and the 
mode- locked fraction as Lml/L. Since only the coherent 
n-row regions contribute to Alg , we can link the value of 
Ala to that of Aja by using Ala = LMLdchAja- In this 
expression Aja can be obtained from the measured pin- 
ning frequency via Eq. (2.8) and (2.6). Using7 = ^oB/pf 
the result for Lml is given by 



A/, 



^ML 



Pf 



fpa 0.67Bdch 



(4.1) 



The different field dependencies of Alg and fp mentioned 
above should thus be attributed to an additional field de- 
pendence of Lml- 

We now evaluate Lml using pf for amorphous NbGe 
films [42,43] and first assume an equilibrium lattice spac- 
ing a = ao(~ 1.075^(f>o/-B) with B = fioH- Figure 
2(c) shows Lml normalized by the channel length L vs. 
field (square symbols). As observed, the coherently mov- 
ing fraction is maximum at B ~ 70 mT. This provides 
a clear definition of the matching field Bm for n = 3. 
At Bm, the longitudinal spacing a in the channel should 
obey a = um = oo and the row spacing 6m = ^o/Bmo-m 
is commensurate with the effective channel width, i.e. 
36m = w- However, away from the matching field the 
array will be frustrated (stretched or compressed) due 
to the confinement. In particular, for B < Bm the lat- 
tice spacing a > oq, while for B > Bm, a < ao- The 
maximum possible difference between a and qq would 
be achieved when the row spacing 6 would not change 
with mismatch, i.e. 6(5) = 6m [44]. This would imply 
a = {Bm / B)aM- Inserting this relation for a in Eq.(4.1), 
the result for Lml, shown by the open circles, is slightly 
modified but shows essentially the same behavior as our 
first analysis: upon increasing the frustration, which we 
define as |1 — (B/Bm)\, the spatial extent of regions with 
3 coherently moving rows shrinks progressively. An ad- 
ditional analysis of the n = 2 ML steps which occur at 
lower fields {B ~ 50mT, where the transition n = 3 ^ 2 
takes place), shows consistently that the spatial extent of 
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the 2 row ML regions increases upon further decreasing 
field. 




f p (MHz) 

FIG. 6. Dynamic ordering frequency /c versus pinning fre- 
quency fp. Dotted line: fc = r/p with r ~ 1 • 10~* s. 

Finally we discuss the behavior of the ordering fre- 
quency fc- As shown in Fig. 5(b), fc{B), denoted by (•), 
exhibits a minimum at the matching field and systemat- 
ically increases with mismatch. Similarly to the decay of 
Lml with increasing frustration, this shows that a larger 
mismatch progressively induces more disorder. We also 
find that, when the field is reduced below i? ~ 50 mT, 
where a 2-row configuration first appears, the ordering 
frequency /"^^ for the emergence of the n = 2 ML effect 
decreases. 

Qualitatively, the behavior of fc is similar to that of 
the pinning frequency fp (or Jc). Both are important 
quantities characterizing the random pinning of a sys- 
tem and are not independent, as follows from a double 
logarithmic plot of fc as a function of /p, shown in Fig. 
6. The data are well described by the relation fc = rf^ 
with T ~ 1 • 10~* s, represented by the dashed line. A 
more detailed fit of the data using a power law relation 
fc oc fp yields an exponent a = 2.1 ± 0.1. 



V. ANALYSIS AND DISCUSSION OF DYNAMIC 
ORDERING AND DEPINNING THRESHOLD 
BEHAVIOR 

A. Comparison with the KV theory 

For a proper discussion of the above results, we first 
shortly describe the phenomenological ordering theory of 
Koshelev and Vinokur (KV) [1]. In their study of a 2D 
vortex system with strong random bulk pinning, they 
found that the shaking action due to motion through the 
pinning potential can be represented by a 'shaking tem- 
perature' Tsh which decreases with velocity as Tsh cm 1/v. 
The dynamic ordering transition occurs when the effec- 
tive temperature T+Tgu is reduced below the cqiiilibrium 
melting temperature T^, i.e. when the velocity exceeds 
Vc oc l/{Tm — T). In later work [4,45] it was shown 



that the shaking temperature refers to (bond) fluctua- 
tions transverse to the velocity and that the associated 
ordering at Vc corresponds to so called transverse freez- 
ing, where inter-chain excursions (permeation modes) are 
strongly suppressed. Within the KV theory we can ex- 
press the ordering frequency fc = Vc/a as: 



fc=VWT^ 



luPl 



^laHchkB{T^-Ty 



(5.1) 



with 7„ the mean squared 2D pinning energy multiplied 
by the area of a pin, pf the flux flow resistivity, a the 
lattice spacing and dch the film thickness. 

In our channel system the ordering can also be de- 
scribed as transverse freezing. We observed this in sim- 
ulations (e.g. for w/Idq ~ 3) as a suppression of the 
inter-chain excursions in parts of the channel at suffi- 
cient velocity [37]. However, different from the 2D sys- 
tem considered by KV, these inter-chain excursions and 
the associated shaking temperature now arise from the 
random interaction with the disordered vortices in the 
CE's and a modification of Eq.(5.1) is required. Think- 
ing in terms of bond fluctuations or a Lindemann cri- 
terion, as in [4], it is clear that it is the short wave- 
length ~ ao disorder component in the potential due 
to quenched vortex displacements d in the CE which 
is relevant for the 'shaking temperature'. This compo- 
nent acts only in a range ~ ao/2 from the first pinned 
row [37], therefore shaking of the outer rows should gov- 
ern the transverse freezing. A rough estimate within 
London theory yields an rms-amplitude of the random 
stress near th e edge ~ EceCee with cee the shear modulus 
and £ce — (a/ (|dp)/ao)/(7r-\/3) representing the random 
strain [16,37] . Taking the longitudinal range of a pin also 
to be ao/2, we replace the parameter 7„ in Eq.(5.1) by 
7ce, resulting in: 



7c 



~ {£ceC&6aobodchf(^f- 



(5.2) 



Further it is important to realize the following: the en- 
ergy scale hgTrn in Eq.(5.1) should be regarded as the en- 
ergy for creation of the dislocation pairs that are required 
for plastic motion i.e. ksTm ksTp ~ CQ^aQdch / {^n) , 
with C66 evaluated at the field and temperature of the 
measurement [46]. For our temperature and fields, this 
energy fcsTp is two orders of magnitude larger than the 
thermal energy ksT which we can therefore neglect in 
Eq.(5.1). Hence, the random shaking (a 1/v) in our case 
essentially represents 'cold working' of the moving struc- 
ture. We also anticipate that the energy ksTp should 
depend on the matching condition since a reduction of 
this energy eventually drives the transition to n = 2 or 
n = 4 rows. Therefore we add a mismatch dependent fac- 
tor Ap-. ksTp = ApCQQaQdch / {^tt) , where Ap is assumed 
to be 1 at matching. Taking into account these changes, 
Eq.(5.1) becomes 



fc 



2Ap^oB 



(5.3) 
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Using the experimental parameters with cge = 
$o-B/(167r^oA^) and A(1.9 K) ~ 1.1 /xm, we obtain the 
value of £ce f rom the value of f,. at matching: Sce — 0.025, 
i.e. -y/ (|dp)/ao ~ 0.13. This is in very reasonable agree- 
ment with our estimate (a/ (|dp)/ao — 0.10) near the 
melting temperature in [16]. Prom Eq.(5.3) and the fc 
data in Fig. 5(b) we can also extract the field dependence 
of Ap characterizing the reduction of the defect creation 
energy. The result is plotted in Fig. 7, showing that close 
to mismatch Ap has decreased by an order of magnitude. 




^0 H (ml) 

FIG. 7. (•) The parameter Ap as determined from fc in 
Fig. 5(b) and Eq.(5.3) using Sce = 0.025. (□) The parameter 

A describing the pinning strength determined directly from 
Eq.(5.4) and the measured critical current density, (o) A de- 
termined from Eq.(5.5) using the field dependence of Ap and 

T = 10"* S. 

Next, we turn to the relation between fc and fp. The 
dc-critical current density can be described phenomeno- 
logically by [25]: 

jc = 2AceG/iBw), (5.4) 

where in our case A varies from A ~ 0.015 at matching 
to A ~ 0.04 — 0.05 at mismatch, see the open squares 
in Fig. 7. Further, when we combine Eqs.(5.1), (5.2), 
Eq.(5.4) with Eq.(2.8) and use 7 = ^oB/pf, we obtain 
the quadratic relation fc = Vc/a = rfp observed experi- 
mentally in Fig. 6 with the timescale r given by: 

r^0.5(£eeM)2J^^. (5.5) 

ApCeepf 

Since Pf,CQ& oc B and experimentally we found that 
r ~ 1 0~^ s, i ndependent of field, this implies that 
A (X ■\/l/(rAp). Using this relation with the field de- 
pendent value of Ap, we obtain A as shown by the open 
circles in Fig. 7. The minima in both data are slightly 
shifted and deviations arc seen for B > B„i, but given 
the approximations made, the overall agreement is still 
reasonable. An important physical implication of the 
relation A cx ^yl/{Ap) is that the increase in pinning 
strength away from matching is directly related to the 
reduction of the defect creation energy. In other words. 



the rise in A reflects an effective softening of the array 
in the channel upon increasing mismatch. Through this 
softening it is able to better adjust to the random CE pin- 
ning potential, very much analogous to the mechanism 
responsible for the peak-effect in ordinary superconduc- 
tors [47]. In our case a more detailed picture of the soft- 
ening mechanism is possible. We already mentioned that 
for B > Bm, the chains are longitudinally compressed 
(a < ao). This will facilitate deviations in the transverse 
direction and reduce the energy to (dynamically) create 
interstitials between rows. For B < Bm, the chains are 
stretched, i.e. a > qq. In this case the energy for a chain 
to accept vortices from a neighboring chain is lowered. 
This in turn facilitates a configurational change in which 
the array can better adopt to the CE potential. An ad- 
ditional mechanism for the rise in A for B < Bm is that, 
due to the mismatch, the outer chains will be pushed to- 
wards the CE, leading to an enhanced influence of the 
CE potential. 

To conclude this section, wc want to mention that the 
relation fc oc fp is not specific to our channels, but can 
be obtained as a general result from the KV theory in 
the strong pinning limit. 



B. dc versus dc-|-rf driven state 

So far we have tacitly assumed that the flow behav- 
ior obtained from our rf-dc measurements simply reflects 
that of the dc-drivcn structure. We now discuss to what 
extent the additional rf-current itself influences the be- 
havior. Recent measurements of the rf-impedance [48], 
which is a sensitive probe of ML at small rf-currcnts, have 
shown that on approaching the dc-driven state, i.e. when 
Irf 0, the voltage broadening SVi^i of the fundamental 
ML step diverges. Since SVi^i oc 6fint, this broadening 
reflects fluctuations in the washboard frequency and, via 
/ = v/a. fluctuations in the velocity and in the longi- 
tudinal lattice spacing a [49]. Correspondingly, the di- 
vergence of SVi^i implies that the dc-driven state lacks 
temporal coherence. Wc also did not observe any nar- 
row band noise in the voltage spectrum of the dc-driven 
state (even for only 30 channels at large dc-drive). At the 
same time, Vi^i oc n remains constant for — > 0, from 
which we conclude that the dc-state still exhibits local 
regions organized in 3 moving chains. Thus, at sufficient 
velocity the dc-state would correspond to temporally in- 
coherent, confined smectic regions [3] of finite (mismatch 
dependent) length, with liquid-like intra chain order and 
residual inter-chain excursions. 

In presence of rf-current the fluctuations are strongly 
reduced, as also observed in experiments on CDW's [50]. 
In simulations we observed that the suppression of the 
inter-chain excursions plays an important role in this pro- 
cess, causing transversely frozen regions in the channel. 
However, the rf-dc IV curves always show incomplete 
ML with the same broadening SVi^i as discussed above. 
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This broadening is too large to be explained by the elastic 
theory in [26]. Therefore, it is either caused by residual 
slip in the n-row regions or by remaining plastic regions 
with interconnecting rows, but further experimental and 
numerical work is required to decide on this issue. 

Finally, we shortly return to the frequency dependence 
of the ML current width /S.Imax in Fig. 4. Extending the 
relation Almaxif) oc Ajmax{f)LML{f) to frequencies 
below fp and taking the ideal tanh(/) behavior Eq.(2.9) 
for Ajmaxj we find that the ordering frequency fc would 
mark the velocity where Lml — > 0. Such interpretation 
implies that the dynamic ordering in our disordered sys- 
tem is a smooth, second order dynamic phase transition. 



VI. SUMMARY 

Using mode-locking experiments, we have investigated 
the dynamics of vortex arrays confined in disordered 
mesoscopic channels. The ML effect allows to trace in 
detail structural transitions from n— l^n— + l con- 
fined moving vortex chains on changing field. A study of 
the amplitude and frequency dependence of the ML steps 
and comparison to simulations of an elastic chain provide 
a complete characterization of the pinning strength, dy- 
namic ordering velocity and coherency of the arrays. We 
find that the spatial extent Lml of coherently moving 
n row regions is large at a matching field and shrinks 
with increasing mismatch. At the same time both the 
pinning frequency fp oc jc and the ordering frequency 
fc (proportional to the ordering velocity) increase with 
mismatch. We show that fc oc fp. Together with our 
previous observation of a divergence of fc near the melt- 
ing temperature in [16], these results provide detailed 
experimental evidence for the phenomenological ordering 
theory of Koshelev and Vinokur [1] . 
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